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Abstract 

For slowly rotating fluids, we establish the existence of a critical point similar 

to the one found for non-rotating systems. As the fluid approaches the critical 

point, the effective inertial mass of any fluid element decreases, vanishing 

at that point and changing of sign beyond it. This result implies that first 

order perturbative method is not always reliable to study dissipative processes 

ocurring before relaxation. Physical consequences that might follow from this 

effect are commented. 
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I. INTRODUCTION 



An alternative path to the study of self-gravitating systems, which avoids the use of 
numerical procedures and/or the introduction of too restrictive simplifying assumptions, 
consists in perturbing the system, compeling it to withdraw from equilibrium state. Then, 
evaluating it after its departure from equilibrium, it is possible to study the tendency of 
the evolution of the object. This is usually done following a first order perturbative method 
which neglects cuadratic and higher terms in the perturbed quantities. This applies when- 
ever the relevant processes occuring in the self-gravitating object take place on time scales 
which are of the order of, or smaller than, hydrostatic time scale. In this case the quasistatic 
approximation fails Jl| (e.g. during the quick collapse phase preceding neutron star forma- 
tion) and the system is evaluated immediately after its departure from equilibrium, where 
immediately means on a time scale of the order of relaxation times. 

Recently, it has been shown that, for systems out of quasi-static approximation, 
a first order perturbative theory is not always satisfactory. In fact, there exist systems for 
which this method seems to be inadequate however small the perturbation is. These ones 
are those for which the parameter 

kT 

a = — - 

r{p + p) 

is close to, or beyond the so called critical point (a = 1). This combination of the tem- 
perature T, the heat conduction coefficient k, the relaxation time r, the energy density p, 
and the pressure p, has been found to be the same in spherically symmetric systems P|§, 
and axially symmetric systems with reflection symmetry ||. Also the viscous spherically 
symmetric case has been studied with similar results ||. 

The astrophysical interest of the study of relativistic rotating fluids is past all doubt. 
Therefore, it seems important to establish, for such systems, the existence, or not, of a 
critical point as described above. With this aim, we assume that, initially, a non-viscous 
slowly rotating object is close to hydrostatic equilibrium (along the r coordinate) and nearly 
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to thermal adjustment (the so called complete equilibrium [p], p. 66]), as measured by a local 
Minkowskian observer. Therefore, the time derivatives of the radial velocity and heat flow 
can be neglected. At that time, we perturb the radial velocity and the heat flow, and 
we evaluate conservation equations and heat transport equation just after the perturbation 
takes place, neglecting cuadratic and higher terms in the perturbed quantities. Here just 
after the perturbation means on a time scale which is of the order of the relaxation time. 
This is necessary if the relevant processes take place on time scales which are of the order, 
or smaller than, hydrostatic time scale. This meaning of just after the perturbation implies 
that physical quantities remain unchanged, but not the time derivatives of the perturbed 
quantities. These ones, are still small, but they cannot be neglected since the system is 
departing from the complete equilibrium. 

As has been mentioned above, it is necessary to use a heat transport equation, together 
with the conservation equations, to find out the existence, or not, of the critical point. 
In order to keep clear of inconsistences, the heat transport equation cannot be the well- 
known Eckart one [[],[7| because it assumes a vanishing relaxation time. Furthermore, this 
theory suffers from two importants drawbacks: Non-causality (the thermal signals propa- 
gate at infinite speed), and unstability (all the predicted equilibrium states are unstable). 
Fortunately, there exist well physically founded thermodynamical theories that avoid these 
problems and that can deal with pre-relaxation processes |8|-[iT|. In this work, we shall 
use the Israel-Stewart heat transport equation. Nevertheless, it is important to emphasize 
that, as in 0-§f, the results found are also valid in the context of the Extended Irreversible 
Thermodynamics |10|JTT[ . 



The paper is organized as follows. The next section is devoted to introduce the interior 
and exterior metrics used, and to construct the stress-energy tensor. Also the validity of 
the slow rotating limit is discussed. In section three, the conservation equations and heat 
transport equation are evaluated just after perturbation and we find the expression for the 
critical point. Finally, we discuss the results in the last section. 

We adopt metric of signature —2 and geometrised units c = G = 1. The quantities 



subscripted with r± denote that they are evaluated at the surface of the object, whereas a 
partial derivative with respect time is denoted by subscript , 0. 



II. ENERGY-MOMENTUM TENSOR 



We consider a nonstatic and axisymmetric distribution of matter and radiation. Let us 
assume that the interior metric is given by [l^ 



ds 2 = Y 2 du 2 + 2—dudr + 2a sin 2 9 (— - Y 2 ^j dud<j> - 2a sin 2 9 ) dnh > 



R 2 d9 2 - sin^ 9 



r 2 + a 2 cos 2 9 + 2a 2 sin 2 9 



Y Y 2 



(1) 



X 2 

where u = x° is a timelike coordinate, r = x l is the null coordinate and 9 = x 2 and <ft = x 3 
are the usual angle coordinates. Here R 2 = r 2 + a 2 cos 2 9, a is the angular momentum 
per unit mass in the weak field limit -the Kerr parameter-, and X and Y are arbitrary 
functions of u, r and 9. The w-coordinate is related to retarded time in a flat space-time 
and therefore, w-constant surfaces are null cones open to the future. In these coordinates 
r-constant surfaces are oblate spheroids. 

The energy-momentum tensor may be expressed in the above coordinates ([!]). Neverthe- 
less, the physical quantities appearing in it will be those measured by a local Minkowskian 
observer comoving with the fluid. Thus, it is necessary to introduce the local Minkowski 
coordinates (t, x, y, z) related to these ones by 



dr / 1 

dt = Ydu + — + a sin 2 9 I — - Y 

dr a sin 2 9 

dx = — H — — d(p 

X X 



(2) 



dy 

dz 



2 , 2 2/1 

r + a cos 9 



1/2 



d9 



2 i 2 2 a] 1 / 2 - 

r + a cos 9 sm 



The radial velocity of matter is given by 

dr XYRujcr — Yau z sin 9 
du 



R(l — uj x ) + Yauj z sin#' 



(3) 



and the orbital velocity is 

n = ^ = ^ , (4) 

du Rsm9 (1 — uj x ) + Yau) z sin 2 Q'' 

where u x and u z are the corresponding components of the velocity of a fluid element as 
measured the locally Minkowski frame. 

In the slow rotating limit a « 1, and consequently Q « 1. Thus, from (|^) 

QRsm8 (1 - uj x ) _ Qr sin6» (1 - u x ) //Vn2 x 
W '~Y(l-na^B)- Y + ° (n) > 



and uj z is also much less than unity. Note that in the static case [i.e. u x = 0) and for a 
local Minkowskian observer (i.e. Y = 1), this means that every fluid element must move at 
non-relativistic velocity [TJ[]. A simple calculus shows that this condition is accomplished 
by most of the known pulsars [[14], p. 146]. 

The interior metric ([!]), can be matched to the Kerr-Vaidya exterior metric [[L^ 

, 2 / 2mr\ , 2 , 4mrasin 2 6 l , , , . 2 , , , 

as =1 — aw + zdudr H ^ a-ua0 — 2a sin Udrdm 

Kb 2 ) r 2 



- R 2 d6 2 - sin 2 6 



9 n 2mra 2 sin 2 ^ 
r 2 + a 2 + - 



(5) 



i? 2 

where R = r 2 -\-a 2 cos 2 ^, a is the exterior Kerr parameter and m is the total mass. It is 
worth mentioning at this point that the metric above is not a pure radiation solution and 
may be interpreted as such only asymptotically |16| . A pure rotating radiation solution may 



be found in reference |[L7| . However, although the interpretation of the Carmeli-Kaye metric 
is not completely clear, the appearance of the critical point is independent of the shape and 
the intensity of the emission pulse, and will be put in evidence for small values of luminosity 
(see below). 



A particular solution can be found in [12|]. Nevertheless, in this work we shall not restrict 
ourselves to a particular solution, and we shall go on using the unknown functions X(u, r, 6) 
and Y(u, r, 9). These ones are constrained by the following junction conditions at the surface 
(r = n) m 



X r=ri Y r=ri I 1 



2mr\ 1/2 



'dX\ fdY\ ( mra 2 sin(2#) \ 



(6) 

/<9Y\ / m (r 4 - a 4 cos 2 + 2r 2 a 2 cos 2 0) \ 

\d^) r=ri ~ \ fr) r=ri = ~ \X{r2- a 2 cos 2 6) (r 2 + a 2 cos 2 #) 2 J r=ri ' 

Next, we assume that, for a local Minkowskian observer, comoving with the fluid, the 

space-time contains: 



mat 



1. An anisotropic fluid of density p mai ; radial pressure p mat and tangential pressure 

2. A radiation field of specific intensity I(x,t;n,v), radiation energy flow q, radiation 
energy density p rad , and radiation pressure p rad . 

The specific intensity of the radiation field I(x, t; n, v), is measured at the position x and 
time t, traveling in the direction n with a frequency v. The moments of I(x, t; n, v) for a 
planar geometry can be written as [TI| 



1 roc rl 
rad 1 - 1 



9 = 2Jo dV J-i dfl I(X} t] V) ' (?) 



1 roc rl 

q — - I dv I dp pl(x,t\n,v) (8) 

2 jo J-i 



and 

I roc rl 



rad = _ 

y 2 



roc rl 

J dv J dp p 2 I(x,t;n, v) . (9) 



where p = cos 8. In classical radiative transfer theory, the specific intensity of the radiation 
field, I(x, t; ft, v) at the position x and time t, traveling in the direction n with a frequency 
v, is defined so that, 

d£ = I(x, t; n, v) dS cos adfl dv dt, (10) 

is the energy crossing a surface element dS, into solid angle dd around n (a is the angle 
between n and the normal to dS), transported by radiation of frequencies (y, v + dv), in 



time dt (see [[y| for details) 



For a nonrotating observer, the radiation portion of the stress-energy tensor reads |T8|JT9 



T. 



R 



p rad _ g 



^rad 



\ 





q p. q 

\(p rad -f ad ) 










\ 



(11) 



/ 



The radiation part of the energy momentum tensor as seen by an observer comoving 



with the fluid can be found by means of a local rotation to (]TT 



T, 



\ 



(12) 



p rad + D 2 p rad _ g q jjq 

-q p rad -Dq 

fl d 

DG -Dq D 2 p rad + p r f d 
where D is an unknown function of u, r and 9 associated with the local dragging of inertial 



V 



/ 



frames effect, fl d = \(p rad - p rad ) and G = \{?>p rad - f ad ). 

The material part of the energy-momentum tensor for this observer is given by 



T ll = (Pm + Pi_)UJJ v - p±r] fll/ + (p- p±)s^s L 



(13) 



where the Minkowski metric is denoted by rj^, s M = 5* and = S t . Thus, the energy- 
momentum tensor, as seen by a Minkowskian observer comoving with the fluid, can be 
written as 



rp rpFi I rpM 



(14) 



In the slow rotation limit D is taken up to first order. Thus, in virtue of (|I2]) and flT3|) , (|I4| 

can be expressed as 



= (p + P±)U fl Uu - P±V^u + {p- p±)spS„ + 1%JJv) + 2q(fiD v ) + 2GU^D 



V), 



(15) 



where q^ = —qs^, = D5^, p = p + p is the total energy density, and p = p + 
prad an( ^ p^ _ p-mat _|_ prad are ^ e total radial pressure and the total tangential pressure 
respectively. 



It remains to express the energy- momentum tensor in curvilinear coordinates ([!]), as seen 
by an observer at rest with respect to the Minkowskian coordinates given by (|2|). Thus, we 
apply a Lorentz boost and the coordinate transformation defined in (0). The boost velocity 
is, in the rotating case, uj = (u x ,0,uj z ) -see for details. Assuming slow rotation limit, 
D, a and uj z are taken up to first order. Thus, 

= (p + p±)U fJ ,U u - p±g fiv + (p~ P±)s^s u + 2q( fl U u ) + 2q^D v) + 2GU^D U ), (16) 



where, is given by ([[]), 



-Y) -cu z r sin9 + 0(ui) 



8t (17) 



+ 



r sin 0^( 7 - 1) + [1 - u x (1 - YX)] + 0{y? z ) 



(18) 



9m = -qsn, 



(19) 



£>„ = 0(uZ)5Z + 0{u 2 z )5l + Dr sin 6 + 0( w 



(20) 



and 



1 
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+ OK 5 



(21) 



Here, 0(u^) corresponds to cuadratic and higher terms in u z , D and a. 



III. DEPARTURE FROM COMPLETE EQUILIBRIUM 



As has been mentioned in the Introduction, we assume that, before perturbation, the 
slowly rotating system is evolving along a sequence of states in which it is close complete 
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equilibrium. Therefore, w-derivatives of oo x and q can be neglected because it is close to 
hydrostatic equilibrium (along the r coordiante) and nearly thermally adjusted. A system 
which is thermally adjusted changes its properties considerabily only within a time scale r c ha 
that is large as compared with the Kelvin-Helmholtz time scale tkh [0 P- 66]. Thus, before 
perturbation we can assume that the w-derivatives of p, p and p± are small, and consequently 
uo x too {i.e. we can neglect cuadratic and higher terms in uj x ). On the other hand, if the 
system is close to hydrostatic equilibrium, then the hydrostatic time scale Th V d ~ \Jr z /m is 
much shorter than the Kelvin-Helmholtz time scale tkh ~ m 2 /2rl, and inertial terms in the 
equation of motion T/ 1 ^ = can be ignored. This condition will be accomplished for small 
values of luminosity I, and consequently for small values of q. Thus, before perturbation 
q~ %). 

We shall evaluate the system just after perturbation. Where, as stated before, just after 
perturbation means on a time scale of the order of the relaxation time. Physically, this implies 
that the perturbed quantities (u x and q) are still much less than unity. Nevertheless, the 
system is departing form complete equilibrium and the w-derivatives of u x and q must be 
small but different from zero (i.e. q t o ~ uj Xj q ~ O(uo x )). 

Thus, the system is characterized by: 

• Before perturbation 

P,o ~ P,o ~ P±,o ~ u x ~ q ~ 0{oj x ) (22) 
~ 9,o ~ 0(u 2 x ). (23) 

• After perturbation 

P,o ~ P,o ~ P±,o ~ u x ~ q ~ uj xfi ~ q )0 ~ 0(u x ) (24) 

In order to clarify the existence of a critical point in slowly rotating fluids, we shall use 
conservation equations (T^ = 0). 
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A. Conservation equations 



Before perturbation, conservation equations read 

R v := = 0. (25) 

After perturbation, physical quantities contained in (|T^) remain unchanged since we are 
evaluating the system on a time scale of the order of the relaxation time. Therefore, the 
only new terms appearing in conservation equations are those containing w-derivatives of lo x 
and q, and conservation equations can be written as 

= R U + u> xfi F v + q fi Qu = 0, (26) 

where tilde denotes that the quantity is evaluated after perturbation, and T v and Q v do 
not depend onuj., q or ^-derivatives of physical variables since we are applying first order 
perturbation theory. The only terms that can contain uj X) q and q$ in T£ = are of the 
form T° . By means of ([[]), (|17|-p0l) and ([16]), we find four equations of the form (^) -see 
appendix A for details- 

R u = (p + p) Wx,o + q,o (27) 
2 

R r = -== Kp + p) w*,o + q,o) 

Re = 
R* = 0. 

Note that R^ is not the total meridional force acting on a given fluid element since it contains 
terms in uj z> q. Nevertheless, i? r does not contain ^-derivatives of physical variables. Thus, 
R r > is the total outward force (pressure gradient + gravitational) along the r-coordinate 
acting on a given fluid element after perturbation. 

The w-derivative of the heat flow q$ can be connected with cD^o by means of an adequate 
heat transport equation. 
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B. Heat transport equation 



As it is well-known, Eckart-Landau transport equation gj^] assumes a vanishing relax- 
ation time. This fact leads to undesirable predictions: An infinite speed for the propagation 
of the thermal signals and unstable equilibrium states [|2T|. Thus, it is necessary to adopt 
a relativistic thermodynamic theory leading to a hyperbolic equation for the propagation of 
thermal signals. On the other hand, we are evaluating the system just after its departure 
form hydrostatic equilibrium and thermal adjustment (in the sense described above). Thus, 
to be consistent with this choice we must use a heat transport equation with non vanishing 
relaxation time. 

We shall use the Israel-Stewart relativistic transport equation [§,[| . For viscous free fluid 



distributions, this one can be written as 22 



rh^U a q u . a + q» = kW v (T u - TU a U u , a ) - ^kT 2 (^^) ^ <f + Tu^q v , (28) 

where k, t and T denote thermal conductivity, thermal relaxation time and temperature 
respectively, = U^U U — g^ u is the projector onto the hypersurface ortogonal to and 
= h^h^U[ a .j3] is the vorticity. 

Before perturbation, transport equations (^) can be symbolized as 

W = 0. (29) 

Just after perturbation (|28"D can be written as 

W + uj x ,^ + q >0 J>* = 0. (30) 

Nevertheless, physical quantities contained on TC 1 do not change just after perturbation, 
and JV 1 = 7"P. Thus, from (^), expresion (|30| ) takes the form 



u x ^ + q fi J^ = Q. (31) 

Vectors X M and as T v and Q u in the preceding section, do not depend on u> x , q or 
^-derivatives of physical variables. 
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The components of ( p8|) containig w-derivatives of q and u x up to first order are 

rh^U a q v , a (32) 

and 

- KTh^U a U v , a . (33) 

Therefore, heat transport equation (|28|) just after perturbation is given by - see appendix B 
for details - 

kT 

q,o = &x,o, (34) 

T 

for any value of 

C. Equation of motion 

We are now in position to find the equation of motion just after perturbation. From (p7[) 
and ( |34| ) we can write 



where 



it = 2( xy P) (l-^)^,o, (35) 



a = —r——:- (36) 
r(p + p) 



As it has been noted in section [III A| , R r > is the total outward force along the r-coordinate 
acting on a given fluid element. Note that it vanishes for a = 1 (the critical point). This 
fact has an important consequence: For a = 1, R r vanishes even though the w-derivative of 
the radial velocity is different from zero. This method also predicts an anomalous behaviour 
beyond the critical point -equation fl35|). If a > 1, then an outward force (R r > 0) implies 
an inward acceleration (cD^o < 0). Therefore, it seems that, for systems out of quasi-static 
approximation, first order perturbation theory can not be applied close to the critical point 
or beyond it. 
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On the other hand, ( p5[) may be compared with the Newtonian form 

Force = mass x acceleration, (37) 

where here, the term 

2 (P + P) n \ /OCX 

XY I 1 - «) (38) 

stands for the effective inertia! mass. Below the critical point, this inertial mass decreases as 
a grows up. This seems to be connected with the dynamical stability of the system, leading 



to a minimum stability for a — 1 [23 



IV. CONCLUSIONS 

In this paper we have studied the departure, of a slowly rotating fluid distribution, 
from a state close to hydrostatic equilibrium (along the r coordinate) and nearly thermally 
adjusted. Our aim has been to elucidate the existence of a critical point similar to the 
found for non- rotating systems [0-0. The existence of this critical point implies that first 
order perturbative method is not always satisfactory to study pre-relaxation processes (i.e. 
processes that take place on time scales smaller than the hydrostatic time scale). 

We have found that, also in this case, there exists such critical point. This one is given 
by condition 

kT 

a = = 1, 39 

T(p + p) 

and it coincides with this one found in spherically simmetric case |2|-^| and in axially simmet- 
ric case 0. Therefore, condition a — 1, establishes an upper limit for which pre-relaxation 
processes can be studied by means of a first order perturbative method. This result is also 
valid if the initial system configuration is strictly in complete equilibrium and radially static 
(i.e. u x = uj xfi = q = q,o = 0). 

Note that this method predicts, for values of a less than unity, that the effective inertial 
mass decreases as a grows. Intuitively, this means that the departure from equilibrium 
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or quasi-equilibrium will be steeper for larger en's, or, in other words, that the smaller a, 
the larger dynamical stability. This point, has been recently illustrated, by means of an 
expression for the active gravitational mass in terms of a |[23|| . It is interesting to emphasize 
that, at least in non-rotating configurations, causality and stability conditions |21[ not always 
forbid the existence of the critical point . 

Finally it is also worth noticing that the critical point and the inflationary equation of 
state for non-dissipative systems (p = —p) are similar in that they imply the vanishing of 
the inertial mass term. Therefore one might wonder about the plausibility of an inflationary 
scenario in a Universe at, or close to, the critical point. 
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APPENDIX A: CONSERVATION EQUATIONS JUST AFTER PERTURBATION 



In the slow rotating limit, U^, s M and read 
7 (!-<*>*) 



Y 



u x X + 0(lu 2 z 



r sin 9 



K (Ai) 



7 (1 - u x ) 
Y 



+ 0{u 
7-1 



r sin 9 \ lu 



8£ + ^X + O^l 



6? 



(A2) 



D 



r sin 9 



+ 0K 2 



6t 



(A3) 



and = —qs^. From fll6|) the only terms that contain w-derivatives of u x and q (up to first 



order in uj x and uj z ) in conservation equations = are of the form T° . In particular 
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(p + P±) (U°U U ) ={ P + PA_) 
(P~P±) (s°s u ) fi = (p-p±) 



7(1 -Ux) 

Y 

7 (1 - u x ) 



U, 



(A4) 



Y 



G (U°D U + U U D°) = G 



D,. 



(q°U u + q°D u + q u U° + q,D°) q = -q fi 



7(1 -Ux) 

Y 

7 (1 - ky) 
y 

7(1- Was) 



(U v + D U + 8 V ) 



Y 



(U u + D U + s u ) 



,0 



Thus, from (|i7| - f20| ), fl2~T|) , (|A4j) and following the definition of T v and Q v given in (p6|) , we 
find up to first order in uo x 



Uxfifu = (p + P_l) 



1 



1 + UJ X 

-w, (p + p) , 



- (p - p±) 



U x fiT r 
^xfiF 9 = 0, 

u x ,oF i> — 0, 





"i - 




V xy ; 


.1 + UJ X _ 


,0 



1 + LU X 

2(p + p) 

xy 



-9 



1 - 

1+UJ X 



(A5) 



q,oQ u = -q,o, (A6) 
r 2 

<?,o£<? = 0, 

where we have neglected in u) x ^T^ terms of the form u Xt0 D, u> Xj oa and u x ^uj z , and in q t oQ^ 
terms of the form q^D, q^a and q t o(d z . Thus, 

(A7) 



.7> = -(p + P) 1, 2/xy, 0, 



and 



Gm = - ( 1, 2/xy 0, 
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(A8) 



The expression of conservation equations can be found by means of (|26|): 

Ru = [(p + p) S^o + q,o] ( 1, 2/XY, 0, ) • (A9) 



APPENDIX B: HEAT TRANSPORT EQUATION JUST AFTER 

PERTURBATION 



The right-hand terms in (|28|) 



2 \kT 2 



q 

;0 



and 



contain factors of the form u X) oq, which are of second order. Therefore, the only terms in 
(28D that contain ^-derivatives of uo x and q up to first order are of the form 



Th^U°q Vt0 (Bl) 



and 



- KTh^U°U u>0 . (B2) 



Using (TAD , (0) and flnjD in flB|) and flB|) 



rh^U q ufi = lil Y UJx) rh^ (-q >0 s v - g Sl/>0 ) (B3) 

Vectors X M and J^ 1 (|3lD do not depend on ct^, g and their w-derivatives because of we are 
using first order perturbation theory. Thus, from (|B3|) 

kT 

u xfl l» = u xfi —W\ (B4) 

and 
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(B5) 



where we have neglected terms of the form q a. Therefore, from ( |B4| ) and (p5|), expresion 



(pl|) takes the form 



9,o — w x q, 

T 



(B6) 



which is valid for any \i. 
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